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ABSTRACT 

In  this  paper  we  axe  interested  in  studying  multiple  de- 
cision procedure  for  k(k  2 2)  populations  which  are  themselves 
unknown  but  which  one  assumed  to  belong  to  a restricted  f£unily. 

We  propose  to  study  a selection  procedure  for  distributions 
associated  with  these  populations  which  eure  convex-ordered  with  ! 

respect  to  a specified  distribution  G assuming  there  exists 
a best  one.  The  procedure  described  here  is  based  on  a statistic  | 

I 

r ' , 

T.  = 7 a.  X.  . for  i = l,...,k  where  X.  . is  the  j-th 

i;].n 

order  statistic  from  F.,  r is  a fixed  positive  integer 

(1  s r s n).  aj  » 9 - 9 g)  for  J = 1 r-1  . : 

a^  = g ^ density  of  G . This  statistic 

T.  was  considered  by  Barlow  and  Doksum  (1972).  If  G(x)  = 

1 

1 - e ^ for  X > 0 , then  nT.  =X.  , +...+  X.  , + 

X i;l»n  i;r-l,n 

(n  - r + 1)  X.  is  the  total  life  statistic  until  r-th 

i;r,n  i 

I 

failure  from  . This  shows  that  the  procedure  based  on  T^ 
generalizes  Patel's  result  (1976)  for  the  IFR  family. 

The  infimum  of  the  probadsility  of  a correct  selection  is 

1 


obtained  and  the  asymptotic  expression  is  also  obtained  using 


the  subset  selection  approach.  Some  other  properties  of  this 
procedure  are  discussed.  We  also  study  the  asymptotic  relative 
efficiencies  of  this  rule  with  respect  to  some  selection  proce- 
dures proposed  by  Barlow  and  Gupta  (1969)  for  the  star-shaped 
ordered  distributions,  Gupta  (1963)  for  the  gamma  populations 
with  unknown  shape  parameters  and  etc.  An  example  is  given  to 
illustrate  the  use  of  the  selection  procedure  for  the  two  inde- 
pendent uniform  distributions.  Application  to  quantile  selection 
rules  for  distributions  convex  ordered  with  respect  to  Weibull 
distribution  is  given.  A selection  procedure  for  selecting  the 
best  population  using  the  indifference  zone  approach  is  also 
studied. 


I 


SUBSET  SELECTION  PROCEDURES  FOR  RESTRICTED 
FAMILIES  OF  PROBABILITY  DISTRIBUTIONS 


1 . I ntroduc  t ion 

In  many  problems,  especially  those  in  reliability  theory, 
one  is  interested  in  using  a model  for  life  length  distribution 
which  is  not  completely  specified  but  belongs,  for  example,  to 
a family  of  distributions  having  increasing  failure  rate  (IFR) , 
or  increasing  failure  rate  on  the  average  (IFRA) . Such 
distributions  form  special  cases  of  what  are  now  commonly 
known  as  restricted  families  of  probability  distributions. 

The  idea  of  using  such  families  stems  from  the  fact  that  in 
many  cases  the  experimenter  cannot  specify  the  model 
(distribution)  exactly  but  is  able  to  say  whether  it  comes 
from  a family  of  distributions  such  as  IFR,  IFRA.  Families 
of  probability  distributions  of  these  types  have  been 
studied  by  several  authors,  see,  for  example,  Barlow,  Marshall 
and  Proschan  [4],  Barlow  and  Proschan  [5,6]  and  Barlow  and 
Doksum  [ 1]  . 

In  this  paper  we  are  interested  in  studying  multiple 
decision  procedures  for  k (k  > 2)  populations  which  are 
themselves  unknown  but  which  are  assumed  to  belong  to  a 
restricted  family.  We  now  give  some  definitions  of  interest 
to  us  (see  Barlow  and  Gupta  [3]). 


(i)  F is  said  to  be  convex  with  respect  to  G 

(written  F -r  G)  if  and  only  if  G ^F(x)  is 
c 

convex  on  the  support  of  F . 

(ii)  F is  said  to  be  star-shaped  with  respect  to  G 

(written  F < G)  if  and  only  if  F(0)  *G(0)  “0 

-1  * 

G ( X ) 

and  ^ is  increasing  in  x i 0 on  the  support 


of  F . 


If  G(x)  = l-  e^,  xzO,  then  F < G is  equivalent 

c 

to  saying  that  F has  increasing  failure  rate  (IFR) . Again 
if  G(x)  =1-6^^,  xsO,  equivalent  to  saying 


N 


that  F has  increasing  failure  rate  on  average  (IFRA) . 

In  the  statistical  literature,  selection  problems  for 
restricted  fzunilies  were  first  investigated  by  Barlow  and 
Gupta  (3].  Some  further  results  in  this  direction  and  a 
review  of  some  important  results  concerning  inequalities 
for  restricted  families  and  problems  of  inference  for  such 
families  have  been  given  by  Gupta  and  Panchapakesan  (10,11] 
and  Patel  ( 15] . 

In  Section  2,  we  propose  and  study  a subset  selection  rule  for 
distributions  which  are  ^ ordered  with  respect  to  a specified 
distribution  G assuming  there  exists  a best  one.  Some 
properties  of  this  rule  are  discussed.  The  infimum  of  the 
probability  of  a correct  selection  is  obtained  and  an 


X. 


■ -yvn-'-  ■ - ■ 

^ • 

asymptotic  expression  is  also  given.  We  also  study  the 
asymptotic  relative  efficiencies  of  this  rule  with  respect 
to  some  selection  procedures.  Section  3 deals  with  selecting 
the  best  population  using  the  indifference  zone  approach.  In 
Section  4,  we  propose  a selection  procedure  for  distributions 
that  eure  < ^ ordered  with  respect  to  G . 

2.  Selection  rules  for  distributions  < ordered  with  respect 

c 

to  a specified  distribution  G . 

Before  discussing  the  selection  problem,  we  give  some 
preliminary  known  results  for  sake  of  completeness.  Let  Q 
be  the  class  of  absolutely  continuous  distribution  functions 
G on  R with  positive  emd  right- (or  left-)  continuous  density 
g on  the  interval  where  0<G<1  and  let  3 be  the  class  of 
distribution  functions  F in  0 such  that  F(0)  = 0 . For 
G €Q  , we  take  g”^(0)  emd  g”^(1)  to  be  equal  to  the  left 
hamd  and  right  hand  endpoints  of  the  support  of  G.  For  F € 5 , 
we  define  f”^(0)  = 0 . For  F€5  and  G€Q  , consider  the 
following  trems format ion  (see  Barlow  and  Doksum  [1]  ) 

(2.1)  = J®"  ^^^g[G"^F(u)]du  , Ostsl  , 

0 

where  g denotes  the  density  of  G . 

We  assume  that  G is  always  fixed.  Since  (the  inverse 

of  Hp ) is  strictly  increasing  on  [0,1]  » ^ is  a 


4 


distribution.  Barlow  and  Doksum  [1]  have  shown  that 


6 if  and  only  if  is  convex  on  the  interval  where 


0<  1 . 


-1 


Since  G is  assumed  known  we  can  estimate  by  substituting 


the  empirical  distribution  of  F ; that  is 


(2.2) 


1 1 1 

“ ! gCG’^F  (u)3du 

no  " 


and 


(2.3)  h;;  (£)  = f'”  9tG- 

is  the  i-th  order  statistic  in  a san^le  of  size 


where  X. 


i,n 

n from  F and  X 


0,n 

(2.3)  can  be  written  as 


0 . If  G(x)  = 1 - 


-X 


for  X 2 0 , then 


(2.4) 


h"^(-)  = ^X,  +...+X  , +(n-r  + l)X  1 

n n n**  l,n  r-l,n  r.n-* 


We  say  that  X,  + . . . + X . + (n  - r + 1)X  is  the  total 

■*  l,n  r-l,n  r,n 


life  statistic  until  r-th  failure  from  F . 


(A)  Selection  procedure  and  its  properties 

Let  , . . . , Ttj^  be  k populations.  The  random  variable 


X^  associated  with  has  distribution  function 


F^,  i ■ 1,  2,  ...»  k , where  f g (i*l»  ...«  k) 


Let 


F.  - denote  the  cumulative  distribution  function  (c.d.f.)  of 
L J 


the  "best"  population.  We  assume  that  (a)  F (x)  sF-  .(x)  for 

1 L ^ J 


J 


all  X , i * 1»  ...»  It  and  (b)  there  exists  a distribution 
G such  that  F.  <G  , i»l,  where  denotes  a 

partial  ordering  relation  on  the  space  of  probability 
distributions.  We  are  given  a seunple  of  size  n from 
each  (i®!#  ....  h)  . Our  goal  is  to  select  a subset 

from  tlie  k populations  so  as  to  include  the  population  with 
. Let  n = {P  = (Fj^....,F^)  :a  a j such  that 
F^(x)  5;Fj(x)  for  all  X and  i = 1,  2,  . . . , k]  . Let 


(2.5)  • 

j=l 


(2.6) 


T = Y a.Y. 

j4i  ^ 


where  X.  . is  the  j-th  order  statistic  from  F.,y.  is 
i;D.n  1 D.n 

the  j-th  order  statistic  from  G,  r is  a fixed  positive 
integer  (1  <r  <n)  , 

a . = go”^  (■^^)  - gc”^  (-^)  for  j = 1, . . . , r-1  and  a = gG  ^ (■^^) 
j n n r n 

For  selecting  a subset  containing  propose  the 

selection  rule  as  follows: 

R^:  Select  population  if  and  only  if 


(2.7) 


T.  2 c,  max  T. 
1 1 • 1 3 

l«3<k 


where  c^^  = c^(k,P*,n,r)  is  the  largest  number  between  0 
and  1 which  is  determined  as  to  satisfy  the  probability 


requirement 


u 


(2.8) 


inf  PtCSlR, ) > 

n ^ 


where  CS  stands  for  a correct  selection,  i.e.,  the  selection 

of  any  subset  which  contains  the  population  with  distribution 

F.,  . . Let  T,,.  be  associated  with  F- . _ and  let  W.  (x) 

[k]  (i)  LiJ  ^ 


be  the  c.d.f.  oi 


(i)  • 


T.gmma  3,1.  Let  F^^,  F^  be  two  distribution  functions  such 
that  F,  (x)  sF-(x)  Vx  and  T.  = Y b X i*l,2,  where 

1 2 i.Zj3i;],n 

b . > 0 for  j c A,  A c f 1 , 2 , . . . , n)  and  X . is  the  j -th 

j ^ 3-  ? J » *• 

order  statistic  from  F^,  i = l,2  , then 

i P[T2  sx]  . 

Proof. 


1 if  T.  > X 
Let  i {o  otherwise 

where  X.,,...,X.  are  n observations  from  F.  (i  = l,2)  . 
il  in  1 

Since  A (X.,,...,X.  ) is  nondecreasing  in  each  of  its 
’ll  in 

arguments,  it  follows  by  induction  (Lehmann  [11]  P.  112^ 
that 


E* 


'*11' 


< E,  (Xjj, 


•'*2n> 


That  is 


pCTj^  s x]  « P[T2  s x]  . 


This  proves  the  lemma. 
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W«  now  mtMtm  and  prove  the  following  theow  which  Is  nore 
general  then  that  of  Patel  £ 15^ . 


Theoren  2.1.  If  €9,  G ^ Q,  P^(x)  fs  and  l-l,2«...,k« 

'fkj  s ® ' *1  * “ *“  1 ' • “■**”  ‘ “ • 

9 8~^(0)  el  and  a^zc^  « then 


(2.9) 


fCCSlRjl  z J 


-1  1 

G (0)  ^ 


where  G_(x)  is  the  c.d.f.  of  T . 
T 


Proof. 


PtCS|Rj]  - PtTj^i  *V(1)  ' •‘■*1 


- I !iA  <f 

0 1 

> I if-)  (*)  (»ir  hmmm,  2.1) 

n 1 


- PCZ^  Z C^z^,  j -1 h-l] 


r -1 

Since  . a.  =gG  (0)<1,  a.  sOvj=l,...,r,  and  (i)(0)  4 0 , 

Li  j ] 

j“l 

by  Lenuna  4.1  of  Barlow  and  Proschan  [51  and  (2.10),  then 


(2.12) 


Since 


^ J y ^ ^(x* . ) . 

1 i;D.n 

1 1 ^ 

~a  2l,  — )a  >l  for  i = 1, . . . , r , 

r C,  j 

i 1 j=l 


and  (p(0)  s 0 


by  Lenuna  4.3  of  Barlow  and  Proschan  [5]  and  (2.10),  we  have 


(2.13) 


(2.14) 


1 ^ 


i 1 j=i 


cp(X^  . ) = Y.  . 

i;],n  St  i;3,n 


where  Y.  . is  the  j-th  order  statistic  from  G , i = 1, 2,  . . . ,1^. 
1 » 3 » n 

Thus  from  (2.11),  (2.13,  (2.13),  and  (2.14), 


P[  max 


l4i<k 


r r 

Z.]  i P[  y a .y,  . j.  c,  y a .Y.  . 

1 L j lc;j,n  1 L j ij],n 


i = 1,2, . . . ,k-l] 

= j 

g’-^(O)  ^ 

This  completes  the  proof. 

Remark  2.1.  (i)  If  g is  nonincreasing,  g{0)  s g G \o)  s 1 and 


“ r * 1 

g G \ — 5 c^,  then  these  conditions  a^j  0,  j = l,...,r,  G (0)  s 0, 
g(r^(0)s  1 arid  Theorem  2.1  are  satisfied. 

(xi)  inf  P[CSIr^]  = r"  ^)dG^(X)  if  g"^0)  = 0 


The  constant  c^^  = c^^  (k, P* , n, r ) satisfying  (2.8)  is  the 


largest  number  between  0 and  1 determined  by 


9 


"g  ^(0)  ‘ * 

We  now  consider  two  specific  distributions  G(x)  , 

**x 

If  G(x)  = 1-e  , xsO  , then  we  have  following  result  which 

slightly  generalizes  the  result  of  Patel  [15]. 


Corollary  2.1.  If  (x)  (x)  Vx  and  i = 1,  . 

Fp,  - < G , G(x)  = l-  e^,x>0  emd 

[Kj  C 

nsmax{r,  } , then 

1 


(2.15)  Inf  pCcsIr,]  = T )dH(x) 

n ^ 6 ‘'i 

2 

where  H(x)  is  the  c.d.f.  of  a random  variable  with 

2r  d.f. 

Proof.  If  G(x)  = 1 - e ^ then  ~ j = l»2,...,r-l 

and  a = — (n  - r + 1)  . 
r n 

Also  ;7-a  il  iff  . 

Cl  r 1-c^ 

2 

By  Theorem  2.1  and  the  fact  that  2nT  is  distributed  as 
with  2r  d.f.,  the  result  follows. 

If  G(x)  = X for  0<  x<  1 , then  we  have  the  following  result 
which  is  a special  case  of  Theorem  2.1  of  Barlow  and  Gupta  [3]. 


in*  PtcslH^3  = „(-J)r  [/r  ’ 1 


0 i=r 


Actually,  the  condition  ^ ® Corollary  2.2  can  be 

relaxed  to  Fn  ^ < G . 

[X]  * 

We  state  emd  prove  the  following  theorem  about  the  asymptotic 

evaluations  of  the  probability  of  a correct  selection  associated 

with  the  rule  in  the  case  where  r is  so  chosen  that 

r^  (n+l)a<r+  1 , 0<a<l  • This  amounts  to  selecting 

populations  with  large  values  of  the  a-quantile  for  a (and  r ) 

r 

as  defined  above.  In  this  case,  — a as  n -♦  » . Note 

n 

that  the  result  holds  for  all  a . 

Theorem  2.2.  If  F.€5,  G€G  for  all  i = 1,  . . . ,k  and 
^ 

(i)  F.(x)  i F ^ (X)  VX  . i = 1 k , Fj^j  X G, 


(ii)  G(x)  has  a differentiable  density  g in  a neighbor- 
hood of  its  a-quantile  t)  , g(ri  ) ^ 0 and  G ^(0)  s 0 , and 

a a 

(iii)  gG  ^ is  uniformly  continuous  on  [0,1],  G ^(x)  is 
convex  and  there  exists  an  f , 0 < ? < 1 , such  that  for  5 s y < 1 


1-y 


is  nondecreasing  in  y , then  as  n -*  ® 


1“C 

(2.17)  PrcSiR.  ] 2 r"  )(-^)^l  d$(x) 

i LC,  c,  a a aa  J 


where  a = 1 - a and  $(x)  is  the  standard  normal  c.d.f. 
Proof . Since  G ^(x)  is  convex,  we  have 


(2.18) 


PCCSjR  T > PfZ  >c  max  Z ] 
^ ' l<j<k  J 


where  Zj^,...,Z^  are  i.i.d.  with  c.d.f.  snd 


is  the  c.d.f.  of  T 


(k)  • 


By  Theorem  (2.2)  of  Barlow  and  Van  Zwet  [7^  and  condition 


(iii) , 


(2.19)  sup|  r 

xaO  b 


gCG  "•F^^,(u)]dul  -»0  a. 


where  F_,  , is  the  empirical  distribution  of  F,,  - . 

Ck]n  ^ Ck] 

Then  we  have  (see  Barlow  and  Doksum  [1])^  for  n large. 


(2.20) 


Z.  2!  Y. 

1 St  i?r,n 


where  Y.  is  the  r-th  order  statistic  from  H„  and 

x,-r.n 

(the  inverse  of  ) is  defined  in  (2.1).  Now 

F F 

m [1^]  -1 

F-.  - "<  G •»  Hp  is  convex.  Since  G (x)  is  increasing  and 

L J c [k] 

convex,  it  follows  that  G H (x)  is  convex.  Since 

.1 

H -<  G and  G (0)  <0  , then  H •<  G . In  a manner 

C^l  c ^[k]  * 

similar  to  the  theorem  (2.1)  of  Barlow  and  Gupta  [3],  we  have 


(2.21)  PCY,  src,max  Y.  2 c,  Y*  , i?^k] 

^ k;r,n  iTr,n-'  k;r,n  • i;r,n  ■' 


where  Y*  is  the  r-th  order  statistic  from  G , i = l,...,k 

i;r,n 

From  (2.18),  (2.20),  (2.21)  and  using  the  fact  that 


(2.22) 


Y*  ~ N(n  , — 

1 ; Y » n Vf  2 , . 

ng  (n  ) 


12. 


the  theorem  follows. 

Before  we  discuss  some  properties  of  the  selection  rule 

, we  introduce  some  definitions  (See  Santner  [16]).  For  a 
given  a(0  < a < 1)  , we  assume  each  has  a unique  a-quantile. 

Let  Fr.T(x)  = F^ . , denote  the  cumulative  distribution  function 
[iJ  [i] 

of  the  population  with  i-th  smallest  a-quantile.  Define 

P„(i)  = is  selected  |r]  where  n...  is  associated  with 

F V^/ 

Definition  2.1. 

(i)  A rule  R is  strongly  monotone  in  if 

jt in  ^[i]  when  all  other  components  of  F^  are  fixed 


P (i)  is 
F 


||  in  when  all  other  components  of  F are  fixed. 


That  means,  P„  (i)  2 P„*(i)  when  F^  . , 2^  F* . , and  P„  (j 

Li  F*  [iJ  St  [1]  F^-’ 


) ^ 


Fr  . T 2 F*  T for 
[ : ] sb  [ ] ] 

j 5^  i . 

Where  F^^  = (Fj.  ^ , . . . , F^  ^ ^ , 

...,F|.j^]),  F* 

IZ  ^ P 

p* 

... F_ . - ) 

p = ^ p 

[1]"' 

' [k]' 

' -2  [1] 

^[j] ^[Jc]^ 

F* 

• ' [j] 

^[k]^  • 

(ii)  A rule  R is  monotone  means 

P (i)  P_(j)  for  all  F^fl  with  F^  . . (x)  2F-.t(x)  . 

£.  F LiJ  LlJ 

(iii)  A rule  R is  unbiased  if  P_(i)  s P„(h)  for  all 

F F 

F Cn  with  Fj.^^(x)  2 F|.j^^(x)  . 


(iv)  A rule  R is  consistent  with  respect  to  n*  meeuis 
inf  PfCSlR]  1 as  n . 

n' 

Theorem  2.3.  If  ^ 0 for  i = l,...,r  , then  Rj^  is 
strongly  monotone  in  . 


Proof. 


Let  ill  (X)  = 


1 if  T...  a c max  T,.. 

i<j<k 


0 otherwise 


where  X 


It  is  similar  to  the  proof  in  Lemma  2.1,  we  can  show  that  R^^ 
is  strongly  monotone  in  rrCi)  . 

Remark  2.2. 

(1)  If  a rule  R is  strongly  monotone  in  tt,..  for  all 

(i) 

i = l,...,k  , then  R is  monotone  and  inf  P[CS1r]  = Inf  P[CS|r] 

" "o 

where  = [F  = (F^,  . . • »F^)  G n ; F^^  = . . . = F^^}  . 

(2)  If  R is  monotone,  then  it  is  unbiased. 

(3)  If  F (x)  = F(x,0.)  , i = l,...,k  and  T.  is  a con- 

1 1 

sistent  estimator  of  9^  , then  R^^  is  consistent  with  respect 
to  n = {F  = (F^ , . . . ,F,  ) ; 3 a j such  that  F . (x)  a F . (x) 

i K 1 j 

for  all  X and  i = l,...,k}  . 

(4)  IfF.€5,G€Q,  F,  <G  , i=  l,...,k  and  the  condition 

1 1C 


(iii)  of  Theorem  2.2  is  satisfied,  we  cam  show  that 


1 


14 


L 


is  consistent. 


The  selection  of  the  population  with  largest  F. 


(i  = l*...fk)  can  be  handled  analogously.  We  assume 


Fj^i^Cx)  <F|.j^^(x)  , i = l,...,k  , and  ® • The  rule 


for  selecting  the  population  with  F.  is  R : Select 

L 1 j z 

population  if  and  only  if 


(2.23) 


c_T.  < min  T. 
^ " lijik  3 


where  c^(0<  c^  « 1)  is  determined  so  as  to  satisfy  the 
basic  requirement.  In  a manner  similar  to  the  proof  of 
Theorem  2.1,  we  have 


Theorem  2.4.  If  F^,  G <:  5 , i]  ^ 1] 


i = 1, . . . ,k 
j = 1, . . . ,r 


F-  , (0)  =0  and  F.  < G and  if  a . i 0 for 


G ^(0)  <0  , gG  ^(0)  si  and  a i c.  , then 

r 2 


(2.24) 


PCcSiR  li  J g’^-^c  x)dG  (X) 

G ^(0)  ^ ^ 


where 


G^(x)  = 1 -G^(x) 


(B)  Efficiency  of  procedure  R^^  under  slippage  configuration. 

Under  the  same  notations  and  conditions  of  Theorem  2.2  emd 
the  comments  above  the  Theorem  2.2,  we  consider  slippage  con- 
figuration F|.^^(x)  = F^'^^  , i = l,2,...,k-l  , and 


Fj.^^(x)  = F(x)  . 0 < 6 < 1 . Let 


xo . 


where  f is  the  density  of  F with  unique  a-quantile,  f 

a 


n (,) 

(2.31)  ARE(R,,R,;6)  = lim  — 

^3 


2 2 

» f {<»  ) 

-0 a_ 

2 2 

a h^(a  ) 
a a 


“X  “X/ 6 

If  G (x)  = l-  e , xvO  and  ~ ^ ® 

2*1  ~ have. 


2.  2 


(2.32) 


ARE(Rj^,R^;ft)  = 


= 0.4803  , n=h  - 


Comparison  with  Gupta  Procedure 


Gupta  [8]  gave  a selection  procedure  for  gamma  populations 


n.'s  with  densities 
1 

The  procedure  R.  is 


x>0,  9.  >0,  i = l,2 k 

r(a)0^  ® ^ 


R^:  Select  population  if  and  only  if 


(2.33) 


X . c , max  X . 

^ 1 • V 3 

l<jr;k 


where  X^  is  the  sample  mean  of  size  n from  and 

c.  is  the  largest  constant  (0<  c.  < 1)  chosen  so  that 
4 4 

P[CS|R.]  . 


For  k = 2,  1]  ~ ^ (see  Barlow  and  Gupta  [3]), 


we  have 


iO. 


n (.) 

(2.34)  ARE(R,,R.:6)  = lim  — — 

^ ^ .^0  "r 


It  is  easy  to  show  that 


2 - 2 
^ a (log  a)  »„  (l-t-6  ) 

2(1-6)^C!;  f(- 
‘a  M 


(2.35) 


ARE(R,,R^;6)  = ARE  ( R,  , R,  ; 6 ) ARE  ( R,  , R^  ; 6 ) 
14  13  3 4 


= f : a log  » Jgi^  Jl  -*•  e \ ^ 

I ./2  (l-6)a  h(a  ) J * 


a 


-X 


If  G(x)=l-s  for  x>0  and  a = l, 
(2.36) 


ARE(R,,R,,M  = _ 


1*4 


(2.37) 


2(l-«)  n 


_ ^~a 


ARE(R, ,R^ ;6t 1)  = 

14  a 


Comparison  of  R^^  and  R^  from  uniform  distribution 

Suppose  and  tt^  arc  two  independent  uniform  populations 

with  distribution  functions  l'\  (i  = l,2)  . 


(2.38) 


x < 0 


F.(x)  =<  — Ocx-'O. 

1 If'.  1 

1 


X > 0 . 

1 


where  6=9 


[11  [2] 


= 1 


A sample  of  n independent  observations  is  drawn  from  each 
of  the  two  populations.  Let  be  the  total  life  statistic 


I 


X > • 


until  r-th  failuro  from  (i=l,2)  where  r < (n+1  )(i " r+1  . 


The  procedure  R is  given  by 


R^;  Select  population  if  and  only  if 


(2.39) 


T*  2 max  T* 

l<3<k 


where  is  chosen  so  that  P[CSlR^‘]  > P*  . Let 


associated  with  i-j  * 


(2.40)  E(SlR^)  -PfCSlR^l  =P[T*^)  '‘''5'^*2)^  *^2^ 


where  Tj^,  are  two  independent  total  life  statistic  until 


r-th  failure  from  uniform  distribution  over  (0,1)  . By  Gupta 

and  Sobel  [ 12] , 


(2.41) 


T ! - u 
1 


N(0,1) 


as  n , 


where  u = 


(2n  - rr  n + 1)  _ ^ np.  (2-gJ.  ^2  ^ 

2n  + 1 


ng  -fln  T ^ u ’ = , a"  = An  and 


. _ g (1-n  ) (2-a)^  11_ 

^ ~ 4 12 


Hence  - ~ = B ^n  where  B = ^ ^ . From  (2.40),  we 

a ~ 'T  A 


have 


T'  - u’  c T'  - u'  c 

EisjR^  - prcsiRj]  - 


20. 


where  are  i.i.d.  with  N(0,1) 


Hence 


E(SiR  ) -PCCSlR  ] = r“  X - (1  )Byh]dt  (x) 

-(1  -^)aAr 

’ c 

5 

Let  E(S|R  ) - P[CSlR  ] = r>0  , we  obtain 


(2.42) 


_ 1)  /H  = 

'c,  6^'"  72  2 B 


6 c. 


Note  that 


inf  P(CS1r  ] = P[T'  > c T' ■]  , where  T'  and  T'  are  defined  as  above, 
o 3 13^  i Z 


T'  - u' 

p[T‘  >ct;:  = p[-^^ > 

X ^ n 


T'  - u" 

=5  1-^1+ '=5 


- P[Z^  ^ c^Z.,  + (c^  - DEk/n]  where  Z^,  Z^  arc  i.i.d.  with  N(0,1) 


Hence 


-{1  ---  )li/n 
‘'5 


prT;?cTn  = r"  ?(—  x - (1  - — )B/n)d»(x)  =tr 7tx_  _ j 

* * J 


Setting  inf  P[CS|r  1 = <•  , we  obtain 
n ^ 

-(1  )B,/h  = r^(e)  Jl  + ~ 


. - / 2 $ (d*’ 

(1  - c_  )./n  = dVi  + c_  where  D = 

5 5 B 


21 


We  see  that  c. 


1 - and  — ~ 1 . F 

« ./n  Tn 


rom  (2.42) , 


(1  + 


J2  D 

/H 


1.  - [ 1 ♦ /• 
6 


B 


/?r 


n 


.m(i-i  i 


-1 


D 


ill  + n** 

B *•6'^ 


Thus 


(2.43)  n 


, rilLulL 

^ B(l-6) 


1 + 6^  - JT  6»~^(p*)  ^ 2 


From  (2.29)  and  (2.43), 

V 

(2.44)  ARE(R  ,R  ,.)  = Um 

c-*0  R_  a n (a  ) 

5 an 

If  we  assume  that  G(x)  = x for  O''  x<’  1 , then 


(2.45)  ARE(R  R ;6)  = ^ ^ 

® 3 ( 1 -fi ) ( 2 -ft ) ■*‘o 

ARE  (Rj^,  R^ ; % ) is  a decreasing  function  of  r,  and  for  n = , 
it  is  equal  to  0.931. 

Note  that  in  (2.45),  R^  is  based  on  r-th  order  statistic  and 
R^  is  based  on  the  total  life  statistic  with  r-th  failure. 

(C)  Selection  procedure  for  distribution  ^ ordered  with  respect 
to  Weibull  distribution 

Assume  that  the  specified  distribution  G(x)  is  given  by 

a 


G(x) 


= ( 


1 - e ^ ^ for  X 2 0 


for  X < 0 


22. 


1 


where  X >0  and  attention  is  restricted  to  « >1  which  is 
assumed  known.  In  this  case,  we  use  as  our  statistic  where 

r-1 

T*  = ) . + (n-r)x"  , i = l k , 

lii.  i;D,n  i;r,n 

j = l 

(as  before,  X.  denote  the  j-th  order  statistic  from 

3 7 j , n 

, i = l,...,k)  . Since  G{x)  is  convex  with  respect  to 

the  exponential  distribution  if  a ;•  1 and  since  the  convex 
ordering  is  transitive,  the  family  of  distributions  which  are 
convex  with  respect  to  Weibull  > 1)  will  have  IFR  distribution. 
Thus  our  interest  here  is  in  a special  subclass  of  IFR 
distributions.  The  rule  for  selecting  the  population  which 
is  associated  with  F,-.  ^ is  as  follows, 

R^:  Select  population  tt^  if  and  only  if 


(2,46) 


T*  > c.  max  T* 

1 ^11  J 

l'']<-k 


where  c^  (0<'  c <rl)  is  determined  so  as  to  satisfy  the  basic 
o o 

probability  requirement. 

Using  the  fact  that  if  F <-  G and  F(0)  =G(0)  = 0 then 

c 

F < G for  n 1 , where  F is  the  c.d.f.  of  , F(x)  is  the 

ri  c 'll  ” 

c.d.f.  of  X,G  is  the  c.d.f.  of  and  G(y)  is  the  c.d.f.  of  Y.  Also, 

n 

G-^F  (X"  ) is  stochastically  equivalent  to  the  i-th  order 

n.  a I/O 


- X 


statistic  from  G*  (x)  = 1-e  , for  x>0  , where  X,  s 

1 , n 

F . 


X 


n,  n 


are  order  statistics  from 


In  a manner  similar  to  the  proof 


of  Theorem  2.1,  one  can  prove  the  following  theorem. 


Theorem  2.5.  If  F^(x)  i = l,...,k  , 


-XX 


, X > 0 , X > 0 and  « (:■  1 ) 


is  known  and  n>max{r. 


r-1 
1 - c. 


] , then 


(2.47) 


inf  P[CSlR^i  = J*  H^'^(^)dH(x) 


where  H(x)  is  the  c.d.f.  of  a 


random  variable  with  2r  d.f. 


(D)  Selection  with  respect  to  the  means  for  Gamma  populations 


Let  be  k populations  with  densities 

X K 

^'"i  a.-l 

f . (x)  = r x^  e ^ . x>0,  B>0*  «.?!,  i = l, 

i r (a  ^ ^ I 


T.et  F.(x)  be  the  distribution  function  of  n.  , i = l,...,k  . 

i 1 

Wc  are  given  a sample  of  size  n from  o'af.’h  i,  . Lc;!  '!’*  I’c 
total  life  statistic  until  r-th  failure  from  . Let 

a <n.,,  be  the  ordered  values  of  « . ' s . We  are 

[11  C’^l  ^ 

interested  in  selecting  the  population  with  the  largest  value 

n ^ 

n . (unknown).  Since  the  mean  of  . is  , selection  of 

[k]  IB 

the  population  with  largest  mean  is  equivalent  to  selecting 


the  population  with  largest  value,  f’pvi  * subset  selection 

L 


rule  based  on  T.  is; 

1 


: Select  population  if  and  only  if 


24. 

(2.48)  T*  > c max  T*  , 

^ l«j<k  ^ 

where  (O--  c^  ' 1)  is  the  largest  value  chosen  to  satisfy 

P[CSl R^]  > P*  . 

Since  the  rule  is  scale  invariant,  we  can  assume 

P = 1 . 

Case  1:  All  are  unknown  and  ? 1 . Let 

n,  = ):  o . ? 1 Yi}  • In  this  case,  by  Corollary 

X i K 1 

2 .2  and  F.  -.^G(x)=l-e  , k > 0 , i = l,...,k  we  have 

^ c 

the  following  result. 

«D 

If  n > max{  r,T~^'l  , then  inf  P[CS|R  ] = f ^(— )dH(x)  , 

7 0^  0 7 i 

2 _ I 

where  H(x)  is  the  c.d.f.  of  a v r.v.  with  2r  d.f.  j 


I 


I 


Case  2 ; n . are  unknown  but  assume  1 < n ■ , i = l,...,k 

1 1 

and  A is  known.  j 

Let  F (x)  be  the  c.d.f.  of  X with  density  function  j 

A 

f (x)  = ■ j'-r  x^'  ^e  ^ , x>0  , 

A T ( A; 

2 

Let  H(x)  be  the  c.d.f.  of  a y r.v.  with  2r  d.f.  and  let  h(x)  be 


(2.49) 


where 


y = F ^{1-e 


P roo  f . 


P[CSlR  ] = PCT*  ^c  max  T*  1 . 


where  T’^ . . is  associated  with  r(_.,  , i = l,...,k  . 

(i)  [i] 

Since  F (x)  <F.(x)  <G(x)  = 1-e  ^ . 

A 1 


(2.50) 


P[CSl 


R_]  > P[T**  ^ c max  T**] 
^ ^ l-ri<k-l  ^ 


where  T**  is  the  total  life  statistic  until  r-th  failure 
k 

from  G(x)  and  T^* ( j = 1 , . . . ,k-l)  is  the  total  life  statistic 

until  r-th  failure  from  F (x)  . 

A 

Since  Asl  then  F < G . Let  ra(x)  =G  F (x) 


(2.51)  P[T**^c^T**.  j = l,....k  - 1]  = F[cp(^  T**)  2 

^7 

m(—  T**)  j =l,....k  - 1] 
n 3 

^7 

By  Lemma  4.1  of  Barlow  and  Proschan  [5"]  with  a,  = . . . = a , = 

^ '1  r-1  n 


(n-r+l)c. 


a = 
r 


, a.  = 0 for  ii  r+1  and  'p(X)  = Y where  X(Y) 
1 st 


is  a r.v.  with  distribution  function  F (G)  respectively,  we  have 

1 ^7  1 

(2.52)  PCfp(-  T**)  >cp(—  T**)  . j = l lc-1]  iPC  ep(-  T**)  > 


2n  ^ ~ ^ 


where  Y^  ( j = 1 , . . . ,k-l)  is  a r.v.  with  x with  2r  d.f. 


I 


From  (2.50),  (2.51)  and  (2.52),  we  have 


P[CS1R^]  > x)dB{x)  , where  B(x)  = P[- T**;  ^ xl 

Since  fo(x)  = -;n(l  - F^(x))  . then  co  ^ (x)  = (1  - e . Thus 

B(x)  =P[T**<ncp^(x)]  =H(2ncp  (x))  =H[2nF^  (1-e  )] 


Now 


dB(x)  2nh  (2nv)  -x  , 

^ = — - ,„v  e where  y = F (1-e  ) 

dx  f.(y)  A 


Hence 


f^(y) 

" x)dB(.)  . ? X)  e-dx 

b ^=7  5 ""t 


This  completes  the  proof. 

Let  S denote  the  size  of  the  selected  subset*  The 
expected  value  of  S when  is  used  is  given  by 


(2.53) 


k 

:(S1R^)  =y  >c  max  T*!  . 

' ' l^j<k 


i=l 


Lot  O ' = 


’ = fa."  (a-j^,»*-»f’^)'  l?n^<A  , i-l,...,'kl  • Poi' 


f.  u ' , since  F^  (x)  <F^(x)  (x)  - 1-e  , then 

E(SlR^)  < -k  P[T**>c  max  T**! 

' ^ 2ij<k 

where  T**  is  the  total  life  statistic  until  r-th  failure 

from  F (x)  and  T**  ( j = 2 , . . . ,k)  is  the  total  life  statistic 
A 3 

until  r-th  failure  from  G(x)  . Hence  E(S|R.^)  < 


k P[T**  > c max  T**  ] . Thus 

^ 2<j<k  ^ 


(2.54) 


sup 


E(S1R,)  = X J dS  (X) 

r\  7 


27. 


where  H(x)  is  the  c.d.f.  of  X r.v.  with  2rd.f.  and  S(x)  is  the  c.d.f. 

of  the  total  life  statistic  until  r-th  failure  from  F (X) . 

L 

Remark  2.3.  (i)  We  can  show  that  the  lower  bound  for  case  2 in 

Theorem  2.6  is  less  than  or  equal  to  the  lower  bound  for  case  1. 

(ii)  Now  we  are  dealing  with  the  problem  in  case  2.  Let 

f”  dn  (x)  = P*  , then  C_  can  be  determined.  If 

0 ^=7  ' ^ 

n a max  fr,  (r-1) / (1-c^ ) 1 , then  we  should  use  the  lower  bound 
for  case  1.  If  r < n < (r-l)/(l-c^)  , then  the  lower  bound  for 
case  1 cannot  be  applied.  In  this  case,  we  can  use  the  lower 
bound  for  case  2. 

(iii)  Sometimes,  the  distribution  function  S(x)  which 
is  defined  above  the  remark  2.3  is  hard  to  compute.  From 

E(sIr^)  5 kP[T**  a T>*  , j = 2,...,k]  where 

T**  is  the  total  life  statistic  until  r-th  failure  from  F and 

L 

Tj*  (j=2,...,k)  is  the  total  life  statistic  until  r-th  failure 
from  G(x)  . Using  the  similar  arguments  in  the  proof  of 
Theorem  2.6,  we  can  get 

X 

E(slR^)  r;  k J*  f‘^  ^1-e  dH(x) 

2 

where  H(x)  is  the  c.d.f.  of  a y r.v.  with  2r  d.f. 

In  this  case,  the  upper  bound  of  E(sIr^)  can  be  computed. 


28 


3.  Selecting  a best  population  - using  indifference  zone  approach. 


Let  be  k populations.  The  random  variable 


associated  with  has  an  absolutely  continuous 


distribution  F.  . We  assume  there  exists  a F-,  _ (x)  such 
1 [k] 


that  F-  . , (x)  > F^ , , (— ) for  all  x , i = l,...,k-l  and 
Cl]  6 


fi  (0<  6 s 1)  is  specified.  Let 
(3.1)  a(6)  ={F  = (F^,  . 

The  correct  selection  is  the  choice  of  any  population  which 


F ) : 3 aj  such  that  F.  (x)  > F . (— )y  i ^ j ] 

k 1 j 0 


is  associated  with  F_,  , . We  propose  the  selection  rule 


; Select  population  if  and  only  if 


(3.2)  T.  = max  T.  whore  T.  is  defined  as  in  Section  1.2 

1 1 1 3 1 


We  want  the  P[  CS  | R ] > P*  , for  all  F^  (t  n , where 

O 


P*  (“<  P*<  1)  is  specified. 


Theorem  3.1.  If  F^€5,  G€Q,i  = l,...,k  , 


where  G^{x)  is  the  c.d.f.  of  T . 


Proof.  P[CS|Rq1  max 


(k) 


l<j«k 


Since  F^^^(6x)  s (x)  , i = l,...,k-l  and  by 


Lemma  2.1,  then 


^lil 


P[CS|Rq1  = > ft  — ^ V j / k]  > PCT^^J  ? 6 T*  V j / k] 


where  T*,...,T*_^  , are  i.i.d.  with  c.d.f.  W^(x)  . 


Using  the  same  argument  as  in  Theorem  2.1,  we  have  our  theorem, 


Remark  3.1. 


inf  PCCS|R8]=  r Gj'”^(^)dG  (x)  if  g"^(0)=0 

Q(ft)  6 


For  given  k,  ft,  P*  and  G(x)  , we  can  possibly  find  the 
values  of  the  pair  {n,r)  , {n  > r)  which  satisfy 


(3.4) 


a > 
r 


« and  T , G!‘'^(f)dG  (x)  a P*  . 

G ^0)  ^ ‘ ^ 


If  G(x)  = X for  0<  xc  1 , we  can  always  find  the 
values  of  the  pair  (n,r)  , (n  r- r)  which  satisfy 

k-1 


(n  - 1 \ or 


n-1 


i=r 


r-1 , , . n-r  ^ * 

X (1  - x)  dx  ^ P 


If  G(x)  = 1-e  ^ for  x^O  , we  can  find  the  smallest 
integer  r , say  r^  , which  satisfies 


J 


is  the  c.d.f.  of  a 


k 


r ^(^)dH(x)  ? P*  where  H(x) 

^0  * 


Y 

r-1 


random  variable  with  2r  d.f.  Since  ~ a >1  iff  n>,  , 

6 r 1-6 


r-1. 


we  can  find  the  minimum  n satisfying  n^maxfr,':; } . 

1—6  ^ 


4.  Selection  procedure  for  distribution  < ordered  with  respect  to  Gj 


Let  t('»0)  be  a given  number.  Let  N^(t)  be  the 


number  of  failures  in  time  t of  the  total  n units  on  life  test 


(without  replacement)  from  which  has  a continuous  distribution 


F. 

1 


i = l,...,k  . We  assume  that 


(4.1) 


(t)  ? 1 , i = 1,  . . . ,k  , 


and  there  exists  a F 


^ (x)  such  that  F (x)  ?F  (x) 

L ^ ,J  ^ L ^ J 


for 


all  X , i = l,...,k  . The  correct  selection  is  the  choice 


of  any  population  which  is  associated  with 


be  the  total  life  statistic  unitl  N.(t)-th  failure  from 

1 


population  rr.  . Let  T,.,  be  associated  with  F,  . , . We 

1 (i)  r i] 


propose  the  rule 


(4.2) 


Rq:  Select  tt.  « T.  > c.  max  T. 

^ 1 1 9 , , D 


l-sj  <:k 


where  Cq(0<c  <1)  is  chosen  so  that  P[CSjR  ] ?•  P*  . 


Theorem  4.1.  If  F.(x)  ?F,,.(x)  Y x , i = l,...,k  and 
, 


F-,  - < G , then 
1^1  * 


(4.3) 


p[cs1r-]  > r )dA- (x) 

^ 5 ^ ""9  ^ 


where  (x)  (A^  (x)  ) is  the  c.d.f.  of  total  life  statistic  until 
n-th  (first)  failure  from  G , respectively. 

Proof.  From  ( 4.1),  we  have 

P[CS|  Rg]  = P[T^^j  > Cg  max  > PCT*  ? CgT*  for  j = l,...,k-l] 

T* 

where  is  the  sample  mean  of  size  n from  F,  . _ , i=l,...,)c-l 

n [1] 

and  — is  the  first  order  statistic  of  size  n from  F.,  . . 
n [k] 

Let 


(4.4) 


T**  X.  for  j=l,...,k-l  . 

3 /-I  3^ 


where 

X X . are 

3l  Dn 

i.i.d.  from  F..  - / 

L J 

j = 1, . . . ,k-l  . 

since 

then  T**  > T*  . 

^ s t ^ 

Hence 

(4.5) 

PCT*  >CgT*  , 

j /k]  > PCT*  >CgT** 

. 3 • 

Let  ^p(x)  = G ^f’-T,,(x)  • 

L K J 

(4.6) 

P[T*  ?CgT**, 

""k 

i ^k]  = PCfp(“)  ? 

"9 

rp(^T**),  j/k] 

Since 

(0  is  star  shaped, 

then 

T*  T* 

,1  Xx  A. 

-0  ( — • — ) > — CP  ( ) . 


(4.7)  'o( — • — ) > — ^(“T'  • 

By  Lemma  3.1  of  Barlow  and  Proschan  [5], 


k 
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(4.8) 


< y n 

n j ^=1  ^ 


. ) . 


y 

NoLc  that  Ti(~)  - Y where  Y,  is  the  first  order  statistic 

n s t 

of  size  n from  G and  cp  (X  ) has  distribution  G , 

^ ' n 

j = l,....k-l  , / =l....,n  . Let  Y.  =7  'PCX.  ) , j = l k-1 

3 r=i  ^ 

From  (4.7)  and  (4.8),  the  right-hand  side  of  (4.6)  is  greater 
than  or  equal  to 


=•  j = i K-iD  = 

- > ’'j  . i = i. 

= “ A^”^(— )dA_(x)  . 

6 ^ S ^ 


This  completes  the  proof. 
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